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Abstract 

The gravitational wave solutions obtained from a perturbation about confor- 
mally flat backgrounds in Einstein gravity are investigated. A perturbation 
theory analysis of the Lesame, Ellis and Dunsby results, based on a covari- 
ant approach, shows that for gravitational waves interacting with irrotational 
dust, the equations are linearization unstable. The gravitational wave equa- 
tions based on the Weyl curvature tensor must be solved by non-perturbative 
methods. The significance of this result for gravitational wave calculations 
and experiments is discussed. 
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I. INTRODUCTION 



Perturbations of spatially homogeneous and isotropic universes have been investigated by 
several authors, beginning with the work of Lifshitz [Ij] , Bonnor []2| , Lifshitz and Khalatnikov 
U, Ehlers @, Hawking §, Sachs and Wolfe §, D'Eath and Moncrief §. More recently, 
a fully gauge invariant formulation of perturbation theory has been presented by Bruni, 
Dunsby, and Ellis [f| . By using the lemma of Stewart and Walker |n| , they derived covariant 



equations from the Bianchi and Ricci identities, which are gauge invariant with respect to 
a zero-order conformally flat background, such as the Friedmann, Robertson, and Walker 
model (FRW). Following this work, Lesame, Ellis, and Dunsby |TT| , |i2"H proved that a tetrad 



frame exists in which the shear tensor and its covariant time derivative are diagonalizable, if 
and only if the divergence of the magnetic Weyl tensor vanishes. Moreover, using the Ricci 
constraints they showed that the magnetic part of the Weyl tensor also vanishes. This result 
is applied to a perturbative analysis of the constraint equations and shows that it leads to a 
linearization instability in the perturbative expansions about a conformally flat spacetime. 
The implications for gravitational wave calculations in GR and for future gravitational wave 
experiments is considered in the following. 

II. FIELD EQUATIONS AND IDENTITIES 

We shall assume Einstein's field equations in the form 

R^iu ~ t;9^ R + A 5V = 8ttGT^ u , (1) 

where A is the cosmological constant, and T^ u is the stress-energy momentum tensor of the 
matter, described by a fluid. The Riemann and Ricci tensors are given by 

^ fj,;[va] 2R [tau'U'oi, Rfj,u Rp ua - (2) 

We have for a perfect fluid: 

T tlv = pu ll u u +ph llu , (3) 



where p is the density, p is the pressure, w M is the velocity of the fluid, = —1, and 

hfj, v = g^ v + u^u u is the projection operator into the hyperplane orthogonal to u M : h^ u u u = 0. 
The gradient of the velocity vector is decomposed as 

Up >v = Uyu, + + -h^Q - u^u u , (4) 

where ii^ = u^ v u v is the acceleration, 6 = is the expansion, = u^-^h^h^ — ^h^O 
is the shear, and oj^ u = U[ a -0\h°h^ is the rotation of the flow lines w M . 

The Weyl curvature tensor can be decomposed into the Ricci tensor R^ v and the Riemann 
tensor R pvaf3 : 

Cp.va.fi = Rpua/3 + 9^Ra\v + 9u[ a R[3\fi + -^R9^[a9/3]u, (5) 

and Cpvap = Cfr v ][ a fl. 

The "electric" and "magnetic" components of the Weyl tensor are defined by 

Ep U = E^ v ) = C Pj avf 3 u a u f3 \ (6) 

and 

= H^ IU ) = -e At7 Q/3 C Q , / 3^ (7 M 7 tt ", (7) 

where is the Levi-Cevita symbol, E^ u and satisfy: E^u" = 0, E^ = and 

H^u u = 0, H/ = 0. 

The Bianchi identities read: 

Rpu[a/3;a] = 0, (8) 

while the Ricci identities are given in @. 

III. PROPAGATION AND CONSTRAINT EQUATIONS 

We shall assume that the vorticity, Up U , and the pressure, p, are negligible. Then, the 
Bianchi identities give Maxwell-type equations |5|,[J: 
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h\E va , p h^ - e^u^Hto = yG/j^p,, (9a) 

h%H ua , p h^ - e^uV^ = 0, (9b) 
+ E^O + h\ p e v)a p p u a H^ p - E a {li a v)oi 

-e^ a e upTS u a u p a 0T E' 7S = -AnGpa^, (9c) 
+ H^d — K 1 ^e u ) a p p u a E 1 l3 ' p — H a ^a l/ ) a 

-e mP(7 e vpT5 u a u p <jP T H° s = 0, (9d) 

where we have used = for p = 0. 

The expansion parameter 9 obeys the Raychaudhuri equation: 

6 + \e 2 + 2a 2 + AnGp = 0, (10) 

where a 2 = \<J pv <J pv . The shear tensor satisfies the evolution equation: 
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<V + cr^aO- ^ - 2#e>> + ^ = 0. (11) 
The Ricci constraint equations obtained from (^) are given by 

/^( V - ^ a / a ;/3 ) = 0, (12) 

and 

= -h\h p u a ( ^ 5 e p)Kr5 u\ (13) 



IV. THE BACKGROUND SPACETIME AND THE CONSTRAINT EQUATIONS 

Since the unperturbed background spacetime is conformally flat, the Weyl tensor van- 
ishes, as do the projected fields E^, H pu , the shear tensor and the vorticity tensor u^. 
Moreover, u M = r. p , where r measures the proper time along the world lines. The spacetime 
is homogeneous and isotropic with 3-surfaces of constant curvature. For the FRW universe, 
the metric is 
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ds 2 = -dr 2 + R 2 da 2 , (14) 

where R = R{t) and da 2 is the line element of space with zero, unit positive or negative 
curvature. 

It has been shown by Hawking || that for p = the rotation satisfies: uj = ujq/R 2 , so 
that the rotation dies away as the universe expands, in accord with our assumption that we 
can neglect uj^. 

Let us define 

We can write equation fl9a|) and (|9b|) as 



(div£) M - E fl p a H /3(T = jGft"^, (15) 

and 

(dwH)p-^ a E^ = 0. (16) 
By considering the proper time derivative of the constraint equation (|T2"|): 

^|(^)--^«(^)-] = 0, (17) 



and using a tetrad frame, Lesame et al., [|Tl|] obtained: 



(TiFi = a 2 H 2 = a 3 H 3 , (1 + A + A># = 0, (18) 

where we used the notation <jj = an (no sum), Hi = Ha (no sum) {i = 1,2,3) and A is a 
constant. Moreover, for the diagonal case: 

<Ti — a, o 2 = Act, 

and 

H\ = XHj H 2 — H . 



Also, the trace-free condition: = gives 

<7g = -(1 + A)<7, H 3 = -(1 + X)H. 

For arbitrary A and a 7^ 0, it follows that H = 0. For 1 + A + A 2 = 0, the values of A are 
complex, but in tetrad form: 

H 2 = -(1 + A + A V = 0, \H\ 2 = (1 + A + X 2 )H 2 = 

so both the shear and the magnetic Weyl tensor are zero, which also follows for the case 
when a = 0, since this implies that the spacetime is FRW for which E = H = 0. This leads 
to the Lesame, Ellis and Dunsby theorem flTW^j: 



For irrotational dust, the divergence of the magnetic Weyl tensor vanishes: 
(divif)^ = for a shear tensor that is diagonalizable in a tetrad 
frame, if and only if the magnetic Weyl tensor vanishes, H pv = 0. 

V. LINEAR GRAVITATIONAL WAVE EQUATION 

We shall now expand the fields E pu , and the tensor a pv in a power series with respect 
to a small parameter e << 1: 

E» v = eE^+e 2 E%+0{e% (19a) 
H, v = eH$ + t 2 E {2 } + 0(e 3 ), (19b) 
^ = ^«+eV( 2 )+0( e 3 ). (19c) 

Let us consider the first-order equations resulting from Eqs. fl9~c|) and (pd|): 

J5« + Ef$e + h\,e v)aPp u a H^» = -AnGpa^J, (20) 
H<$ + H<$6 - h\,e v)aPp u a E^ = 0. (21) 

By multiplying fl2§) by d/dr = u a V a and (gj) by /i^ (p e /3 ) 7i/T u 7 V r , we get § 

E% ~ A 2 E« + liS®0 + E®0 + ^ 2 + ^Gp) + 8vrGaW(lp + \p) = 0, (22) 
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where 

A 2 E, U = (E { ^ T h a sh^h T ,), a h^h s ^ u 

is the Laplacian operator in the hypersurface r = constant. It has the eigenfunction expan- 
sion: 

with m = 0. 

For a non-expanding congruence u a and empty Minkowski spacetime, Eq. fl22|) reduces 
to the wave equation 

= 0. (23) 



VI. LINEARIZATION INSTABILITY 



Equating coefficients of powers of e, we obtain from Eqs. ( pa|) and (|9bD: 

(div E^), = ^-Gh%p ;u , (24a) 

(divtfW^O, (24b) 

(divE^), - Ej&ffW" = Y Gh "^ ( 24c ) 

(divF( 2 )),-S( 1 J CT E( 1 )- = 0, (24d) 



Svr 



(div fiW), - £ >" = -G7i> ;i ,, (24e) 

i=i d 

n-l 

(div - X) E« ffJ BM- = 0, (24f) 

8=1 

and from Eqs. fl9c|) and (pd|): 

^ + 4 X > + ^(^, )a/3 pw a ^ 7 (1)/3;p = -^Gpa$, (25a) 
«J5 + - h\,e v)aPp u a EW'» = 0, (25b) 
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E$ + E$6 + h\ p e u)af3p u a H^ - E^ a {p a^K )a 

-e mP(T e vpT5 u a u p a^ T E^ s = -AvGpaf} , (25c) 
H$ + H$6 - h\ p e u)af3p u a E^ - H^ a ^ v)a 

-£^e^VW ri = 0, (25d) 



n-l 

(n-i) 



i=l 

-e paPa e vpT6 u a u» £ ^ = -^Gpa^j , (25e) 

i=l 

n-l 

re— 1 

-e^e upT5 u a u p £ = o. (25f) 



i=i 



By inserting the infinite power series in e for H pv into the result of the Lesame, Ellis 
and Dunsby theorem: (divi^) M = and H pu = 0, we find that vanishes in every order. 
Then, Eq. ( |24b| ) leads to the result: 

flJS = 0, (26) 

and to first-order we get 

(&ivE^) p = ^Gh\p, u , (27a) 

E^+E^e = -ArcGpa^, (27b) 

curl E$ = h\ p e v)afip u a E^ = 0. (27c) 

In second-order, we obtain 

{dwE^) p = ^Gh\p, V) (28a) 

(divff^-E^W^O, (28b) 
E$! + E$9 + h\,e v)a ^E^ - E^\ )a 

-e pa ^e uprS u a u p a^ T E^ s = -4nGpa$ , (28c) 

HjS + H$9 - h\ p e u)aPp u«E^o = 0. (28d) 
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Thus, neither the first-order nor the second-order of perturbation theory predict the same 
evolution as the rigorous theory. This means that there is a linearization instability and 
the gravitational wave equations can only be solved by non-perturbative methods. It is not 
possible to derive the standard linear wave equation fl23|) in empty spacetime, because of 
the absence of a curli^ 1 ) term in the first-order equations. This confirms the Lesame et 
al., conclusions [fLl]| . However, our analysis reveals that the breakdown of the linearization 



scheme follows inevitably from the Lesame, Ellis and Dunsby theorem, when a power series 
solution of the equations is sought. 

VII. CONSEQUENCES FOR GRAVITATIONAL WAVE EXPERIMENTS 

Considerable work has been devoted to studying the linearization stability of gravita- 
tional wave perturbation theory in Minkowski spacetime P, |1"3| , |T4|| . It demonstrated that 



Cauchy stability can be proved and that the linear approximation of the field equations can 
be expected to be stable for the empty space Einstein equations. However, the problem of 
perturbative stability for cosmological solutions, such as the FRW solution of the Einstein 
field equations, is a more difficult issue to resolve. The full Einstein field equations have 
(under certain conditions) a Cauchy solution which is stable, but potential difficulties arise 
when a perturbative analysis of the equations is performed due to the role of the Bianchi and 
Ricci constraints. This is borne out by the results of Lesame et al., who consistently retain 
the Ricci constraints in their analysis by demanding that they hold in the time evolution of 
the initial data. A central issue is that the Lesame et al., analysis is based on a fully gauge 
invariant (covariant) perturbation theory. As first stressed by Hawking ||, a non-gauge 
invariant analysis based on metric perturbations is not reliable and can yield misleading 
results. He therefore suggested that the perturbation theory be based on an analysis of the 
curvature tensors. 

In the proposed gravitational wave experiments such as the LIGO/ VIRGO or LISA 
projects [113], which attempt to detect gravitational radiation directly, the gravitational 
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waves arrive at the detector after travelling through the universe described by a cosmo- 
logical model such as the standard FRW model. In these models there is always a non- 
vanishing density p i.e., the estimated density of the smoothed out fluid of the universe is 
p ~ 10~ 29 gcm~ 3 . Even though this is a small density it will play an important role in 
the perturbative stability analysis, because the linearization instability of the gravitational 
wave equations in the presence of an irrotational dust is expected to be singular, in the 
sense that the solutions of the gravitational wave equations become singular in the limit 
of a conformally flat FRW solution. Thus, only non-pertubative or non-linear numerical 
solutions of the gravitational wave equations, based on a gauge invariant set of equations 
within a cosmological model scenario, can be trusted to produce reliable physical predictions 
for gravitational wave experiments. 

The proposed strong gravitational wave sources which one hopes can produce detectable 
radiation, such as the coalescence of neutron star and stellar mass black hole binaries in 
distant galaxies or colliding black holes, involve strong gravitational fields with non- vanishing 
matter and shear in the close vicinity of the sources, requiring a non-perturbative solution 
of the wave equations. 

Although the conclusions drawn about the linearization instability of the cosmological 
gravitational wave equations by Lesame et al., and in the present work, is based on irrota- 
tional dust, such an approximation is expected to be very reliable in the present universe, 
since the vorticity and pressure are vanishingly small. Therefore, the linearization instabil- 
ity of the gravitational wave solutions travelling through a medium discussed here, plays an 
important role for realistic gravitational wave calculations and experiments in the present 
universe. 

Perhaps, this kind of linearization instability of the gravitational wave equations, arises 
from a breaking of the conformal symmetry of the cosmological model e.g., the FRW model 
which is conformally invariant with a constant curvature. Further work should be carried out 
to understand the mechanism associated with this symmetry breaking, which is responsible 
for producing the linearization instability. This mechanism may be related to the result 
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obtained by Kundt and Triimper and Szekeres [|L6| , |17 l) that in an exact treatment of GR, 
no Petrov type N irrotational dust solutions exist. 



VIII. CONCLUSIONS 

A perturbative analysis of the Weyl curvature Maxwell-type equations and the associated 
constraints about a conformally flat background, which was chosen to be represented by the 
FRW model, led to a linearization instability when the Lesame, Ellis and Dunsby theorem 
was invoked. This confirmed the conclusions about linearization instability inferred by 



Lesame et al., |Tl|Jl4 . 

It is important to emphasize that this instability is inevitable already at first-order in a 
conventional power series expansion of the Weyl curvature equations in a small parameter 
e, and reflects itself in the higher order equations. The Ricci constraints must be imposed 
at the lowest order, which results in Hj^J vanishing and the breakdown of the power series 
perturbation theory. Thus, the perturbative instability is a generic feature of any consistent 
expansion about an FRW backround, based on the covariant Weyl curvature approach. Per- 
haps, a reinvestigation of the Cauchy perturbative instability of gravitational wave equations, 
expanded about an irrotational dust solution, is needed to fully understand the implications 
of this result. 

Since gravitational wave experiments are conducted for strong gravitational field sources 
with non-vanishing matter density and shear, and since the gravitational waves that the 
experiments hope to detect travel through an irrotational dust in the present universe, then 
the whole issue of the validity of the perturbative calculations of gravitational wave solutions 
should be treated with caution. This may become an important topic of study in numerical 
calculations of gravitational waves in the near-linear regime [|1^ . 
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